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Abstract 

The effective permittivity dyadic of a composite material containing particulate constituent 
materials with one constituent having the ability to display the Pockels effect is computed, 
using an extended version of the strong-permittivity-fluctuation theory which takes account of 
both the distributional statistics of the constituent particles and their sizes. Scattering loss, 
thereby incorporated in the effective electromagnetic response of the homogenized composite 
material, is significantly affected by the application of a low-frequency (dc) electric field. 
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1 Introduction 



Two (or more) materials, each composed of particles which are small compared to all relevant 
wavelengths, may be blended together to create an effectively homogeneous material. By ju- 
diciously selecting the constituent materials, as well their relative proportions, particle shapes, 
orientations and sizes, and distributional statistics, the homogenized composite material (HCM) 
can be made to display desirable magnitudes of effective constitutive parameters [1, 2]. Further- 
more, an HCM can exhibit effective constitutive parameters which are either not exhibited at 
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all by its constituent materials, or at least not exhibited to the same extent by its constituent 
materials [3, 4]. A prime example of such metamaterials is provided by HCMs which support 
electromagnetic planewave propagation with negative phase velocity [5]. 

The focus of this paper is on the electromagnetic constitutive properties of HCMs [1, 6]. If 
one (or more) of the constituent materials exhibits the Pockels effect [7], then a further degree 
of control over the electromagnetic response properties of an HCM may be achieved. That 
is, post-fabrication dynamical control of an HCM's performance may be achieved through the 
application of a low-frequency (dc) electric field. For such an HCM, the potential to engineer its 
electromagnetic response properties (i) at the fabrication stage by selection of the constituent 
materials and particulate geometry, and (ii) at the post-fabrication stage by an applied field, is 
of considerable technological importance in the area of smart materials [8]. 

The opportunities offered by the Pockels effect for tuning the response properties of composite 
materials have recently been highlighted for photonic band-gap engineering [9, 10] and HCMs 
[11]. In particular, the predecessor study exploiting the well-known Bruggeman homogenization 
formalism revealed that the greatest degree of control over the HCM's constitutive parameters 
is achievable when the constituent materials are distributed as oriented and highly aspherical 
particles and have high electro-optic coefficients [11]. However, the Bruggeman formalism may 
not take predict the scattering loss in a composite material adequately enough, as it does not 
take into account positional correlations between the particles. Therefore, in the following 
sections of this paper, we implement a more sophisticated homogenization approach based on 
the strong-permittivity-fluctuation theory (SPFT) which enables us to investigate the effect of 
the dc electric field on the degree of scattering loss in electro-optic HCMs. 

A note about notation: Vectors are in boldface, dyadics are double underlined. The inverse, 
adjoint, determinant and trace of a dyadic M are represented as M~ 1 , M , det ( M), and 

tr (M), respectively. A Cartesian coordinate system with unit vectors Ux,y,z is adopted. The 
identity dyadic is written as /, and the null dyadic as 0. An exp(— iut) time-dependence is 
implicit with i = sj— 1, oj as angular frequency, and t as time. The permittivity and permeability 
of free space are denoted by eo and fio, respectively; and the free-space wavenumber is ko = 

2 Theory 

Let us now apply the SPFT to estimate the effective permittivity dyadic of an HCM arising 
from two particulate constituent materials, one of which exhibits the Pockels effect. Unlike 
conventional approaches to homogenization, as typified by the usual versions of the much-used 
Maxwell Garnett and Bruggeman formalisms [1, 6], the SPFT can take detailed account of the 
distributional statistics of the constituent particles [12, 13, 14]. The extended version of the 
SPFT implemented here also takes into account the sizes of the constituent particles [15]. 
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2.1 Constituent materials 



The two constituent materials are labeled a and b. The particles of both materials are ellipsoidal, 
in general. For simplicity, all constituent particles are taken to have the same shape and 
orientation, as specified by the dyadic 
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£ = (aia 2 «3)" 1/3 ^2 a K^K^K, (1) 

K=l 

with cxk > OVK G [1,3] and the three unit vectors aj.,2,3 being mutually orthogonal. Thus, the 
surface of each constituent particle may be parameterized by 

r e (e, ( p) = v g'i(e,<p), (2) 

where r(9, <j)) is the radial unit vector in the direction specified by the spherical polar coordinates 
9 and <p. The size parameter r\ > is a measure of the average linear dimensions of the particle. It 
is fundamental to the process of homogenization that rj is much smaller than all electromagnetic 
wavelengths [1]. However, 77 need not be vanishingly small but just be electrically small [16, 17]. 

Let V a and Vj, denote the disjoint regions which contain the constituent materials a and b, 
respectively. The constituent particles are randomly distributed. The distributional statistics 
are described in terms of moments of the characteristic functions 



*<(r) 



1, r G V t , 

(£ = a,b). (3) 

0, r V e , 



The volume fraction of constituent material £, namely fg , is given by the first statistical moment 
of &e ; i.e., ( $^(r) ) = fg; furthermore, f a + fb = 1- For the second statistical moment of &e, we 
adopt the physically motivated form [18] 



(^(r) <D £ (r')> 



( (^(r))($,(r')), \UT l ' (r-r')| > L, 

(4) 

, <^(r)), l^" 1 ' (r-r')l < L. 



The correlation length L > herein is required to be much smaller than the electromagnetic 
wavelength(s), but larger than the particle size parameter rj. The particular form of the covari- 
ance function has little influence on SPFT estimates of the HCM constitutive parameters, at 
least for physically plausible covariance functions [19]. 

Next we turn to the electromagnetic constitutive properties of the constituent materials. 
Material a is simply an isotropic dielectric material with permittivity dyadic = J in the 
optical regime. In contrast, material b is more complicated as it displays the Pockels effect. Its 
linear electro-optic properties are expressed through the inverse of its permittivity dyadic in 
the optical regime, which is written as [11] 



+ s 4 (u 2 u 3 + U3U2) + s 5 (111U3 + U3U1) + s 6 (uiu 2 + u 2 ui) >, (5) 
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where 



sj=Y, r JK E K, J €[1,6]. 



(6) 



K=l 

The unit vectors 

ui = — (u x cos 4>b + u. y sin <^,) cos #6 + u 2 sin 0& ' 
u 2 = u x sin 06 — u y cos <pb 

u 3 = ( u :r COS (f>b + U y sin <pb) sin 6*fe + u 2 cos Ob 



> , 9b £ [0, tt] , </>6 £ [0, 2tt] 



(7) 



pertain to the crystallographic structure of the material. In (5) and (6), = ur- • E dc , 
K G [1, 3], are the Cartesian components of the dc electric field; ef\ 3 are the principal relative 
permittivity scalars in the optical regime; and rjx, J € [1,6] and K G [1,3], are the 18 
electro-optic coefficients in the traditional contracted or abbreviated notation for representing 
symmetric second-order tensors [20]. Correct to the first order in the components of the dc 
electric field, which is commonplace in electro-optics [21], we get the linear approximation [22] 
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£ 

K=l 



-K 



(P) 

e K 's K I VL K U K 



- 4 6) 4 &) s 4 (u 2 u 3 + u 3 u 2 ) - efeg' s 5 (uiu 3 + U3U1) - ef'e^' s e (uiu 2 + u 2 Ui) } (8) 



,(b)(b) 



from (5), provided that 



{max 1 (b) 1 1 r max i i ~i . . -, 
jre[i,3] \tk \ J \J6[i,6] < 1 • 



(9) 



The constituent material b can be isotropic, uniaxial, or biaxial, depending on the relative values 



of ef\ and e^' . Furthermore, this material may belong to one of 20 crystallographic classes 
of point group symmetry, in accordance with the relative values of the electro-optic coefficients. 

In order to highlight the degree of electrical control which can be achieved over the HCM's 
permittivity dyadic, we consider scenarios wherein the influence of the Pockels effect is most con- 
spicuous. Therefore, the crystallographic and geometric orientations of the constituent particles 
are taken to be aligned from here onwards; furthermore, E dc is aligned with the major crystal- 
lographic/geometric principal axis. For convenience, the principal crystallographic/geometric 
axes are taken to coincide with the Cartesian basis vectors u Xj y tZ . 
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2.2 Homogenized composite material 



The bilocally approximated SPFT estimate of the permittivity dyadic of the HCM turns out to 
be [15] 

f CM = e [dmg{e% CM ,e% CM ,e? CM )] = | cm - (l + g • (10) 

Herein, e cm is the permittivity dyadic of a homogeneous comparison medium, which is delivered 
by solving the Bruggeman equation [11] 

faX (a) + fbX ib) =g, 
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where 

X® = (| W - g cm ) • [l + R- - i m ) ] ~* , (1 = a, b) (11) 
are polarizability density dyadics. The depolarization dyadic is represented by the sum 

R = R° +R >0 (v)- (12) 

The depolarization contribution arising from vanishingly small particle regions (i.e., i] — > 0) is 
given by [11, 23] 

e? = ir ^-L' " e * " - <i3) 



tr (| cm • 



where 



A = diap 



sin 9 q cos<^ g \ / sin sin0 ? \ 
u x • U ' u x J ' \ u„ • U • u„ J ' \ u z • U • u 



COS # n 



2/ ^_ u 9 



(14) 



while the depolarization contribution arising from particle regions of finite size (i.e., rj > 0) is 
given by [15] 



D>°(V) = f- 



3 />27T 



47T 



3 («+-«-) , . ( f | 



a 



+ilo 2 iiq [ k1- kI) (3 >, 



(15) 



with 



[2g cm -tr(V m ) /] • A — tr (f cm • A) I 
det (| cm ) 



tr e cn ]. - A 

-ad] = 



tr (a) tr (4) 



tr (e cm • A) 



= = ad 3 tr (e cm • A) 



A, 



tr \i c Zj ' A ) ~ tr (lad,- j tr (4) J - 4det (i cm ) tr (4) tr (i cm * 4) 

tr (| c a ™. ) tr (A) - tr • ± A 



2tr(^) tr (| cm • 
The mass operator dyadic in (10) is defined as 



-A, 
(16) 
(17) 

(18) 
(19) 

(20) 



3 Numerical results and conclusion 



The prospects for electrical control over the HCM's permittivity dyadic are explored by means 
of an illustrative example. In keeping with the predecessor study based on the Bruggeman 



5 



homogenization formalism [11], let us choose an example in which the influence of the Pockels 
effect is highly noticeable. Thus, we set the relative permittivity scalar = 1, whereas 
the constituent material b has the constitutive properties of potassium niobate [24]: = 4.72, 
= 5.20, 4 6) = 5.43, r 13 = 34x 1(T 12 m V -1 , r 23 = 6xl(T 12 mV" 1 , r 33 = 63.4x 10" 12 m V -1 , 
r 42 = 450 x 10~ 12 m V -1 , r 5 i = 120 x 10~ 12 m V -1 , and all other tjk = 0. For all calculations, 
the volume fraction is fixed at f a = 0.5; the shape parameters are a\ = 012 = 1 and a 3 = 9 (so 
that the particles are prolate spheroids); the crystallographic angles 0b = 4 > b = 0; E\% = an d 
the range for E$ c is dictated by (9), thereby allowing direct comparison of the SPFT results 
with results for the Bruggeman formalism [11]. The SPFT calculations were carried out for an 
angular frequency 00 = -it x 10 15 rad s^ 1 , corresponding to a free-space wavelength of 600 nm. 

In relation to the predecessor study [11], the effects of two quantities have to be considered: (i) 
the correlation length L and (ii) the particle size parameter r\. Let us begin with L. In Figure 1, 
the three relative permittivity scalars e.^y¥ computed from (10) in the limit 77 — > (by setting 
D_ >0 = 0), are plotted as functions of k$L and E^ c . At L = 0, the variations in ^y,^ 1 as 
E^ c increases are precisely those predicted by the Bruggeman homogenization formalism [11]; 
but the imaginary parts of e^y^ 1 are nonzero for L > 0. The emergence of these nonzero 
imaginary parts is unambiguously attributable to the incorporation of scattering loss in SPFT. 
As L increases, the interactions of larger numbers of constituent particles become correlated, 
thereby leading to an increase in the overall scattering loss in the composite material. This is 
most noticeable in the z direction (which is direction of the major principal axis of U and eW, 
and also the direction of E^ c ), as indicated in Figure 1 by the behavior of e^ CM . At k$L = 0.2, 
the magnitude of the imaginary part of e^ CM increases by approximately 50% as ranges 
from 3 xlO 8 to -3 xlO 8 V m _1 . 

The effect of the size parameter r\ is very similar (but not identical) to that of the correlation 
length, as is apparent from Figure 2 wherein the real and imaginary parts of are plotted 

against k^r] and E$ c , for L = 0. The nonzero imaginary parts of ^y ^ 1 , which arise for 77 > 0, are 
also attributable to coherent scattering losses associated with the finite size of the constituent 
particles [17]. The effect of the dc electric field over the real and imaginary parts of the HCM's 
relative permittivity scalars at a given value of r\ (with L = 0) is much the same as it is at the 
same value of L (with i] = 0). 



Lastly, let us turn to the combined influences of the correlation length and the size of the 

HCM 

x,y,z 



constituent particles. In Figure 3, plots of the real and imaginary parts of e^^^ 1 against korj 



and Ef c are displayed for the case where k$L = 0.25. Plots of the real parts of e^y ^ 1 are similar 
to those in Figures 1 and 2. In contrast, the magnitudes of the imaginary parts of t^ylz 1 are 
substantially larger in Figure 3 than in Figures 1 and 2. As is the case in Figure 2, the imaginary 
part of e^ CM increases by approximately 50% as Ef 2 ranges from 3 xlO 8 to —3 xlO 8 V m _1 , 
at korj = 0.2. 

Let us conclude by stating that a particular feature of the electromagnetic response of an 
HCM is that attenuation can arise, due to coherent scattering loss, regardless of whether the 
constituent materials are dissipative or nondissipative. The extended SPFT [15] provides a 
means of estimating the effect of this scattering loss on the effective constitutive parameters 
of the HCM, and relating it to the distributional statistics of the particles which make up the 
constituent materials and their size. When one (or more) of the constituent materials exhibits 
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the Pockels effect, the degree of scattering loss may be significantly controlled by the application 
of a low-frequency (dc) electric field. The technological implications of this capacity to control 
dynamically, at the post-fabrication stage, the electromagnetic properties of an HCM — which 
itself may be tailored to a high degree at the fabrication stage — are far-reaching, extending 
to applications in telecommunications, sensing, and actuation, for example. 
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Figure 1: Real (left) and imaginary (right) parts of the HCM's relative permittivity scalars 
e x,yf> as estimated using the SPFT, plotted against E$ c (in V uT 1 xlO 8 ) for Ef% = and 
k L. The constituent material a is air and the constituent material b is potassium niobate. The 
particles of both materials a and b are parallel prolate spheroids, characterized by the shape 
parameters a\ = ol-z = 1 and 03 = 9. Calculations were made in the limit 77 — > 0. 
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Figure 2: As Figure 1 but with the real (left) and imaginary (right) parts of the HCM per- 
mittivity scalars ^x,y,¥ plotted against E^ c (in V in -1 xlO 8 ) and kor], with correlation length 
/. 0. 
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Figure 3: As Figure 1 but with the real (left) and imaginary (right) parts of the HCM permit- 
tivity scalars txyjf plotted against E$ c (in V mT l xlO 8 ) and karj, with k L = 0.25. 
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